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Abstract 

In  recent  years,  researchers  have  been  dedicated  to  the  study  of  underactuated  manipulatcrs 
whieh  have  more  joints  than  control  actuators.  In  previous  works,  one  always  assumes  that 
there  is  enough  dynatnie  coupling  between  the  aetive  and  the  passive  joints  of  the  manipula¬ 
tor,  for  it  to  be  possible  to  control  the  position  of  the  passive  joints  via  the  dynamic  cou¬ 
pling.  In  this  work,  the  authors  aim  to  develop  an  index  to  measure  the  dynamic  coupling,  so 
as  to  address  when  control  of  the  underactuated  system  ts  possible,  and  how  the  motion  and 
robot  configuration  can  be  designed.  We  discuss  extensively  the  nature  of  the  dynamic  cou¬ 
pling  and  of  the  proposed  coupling  index,  and  their  applications  in  the  analysis  and  design 
of  underachjated  systems,  and  in  control  and  planning  of  robot  motion  configuration. 


1 


1  Introduction 

In  recent  years,  researchers  have  been  turning  their  attention  to  so  called  underactuated 
systems,  where  the  term  mderactuated  refers  to  the  fact  that  the  system  has  more  joints  than 
control  actuators.  Some  example.';  of  underactuated  systems  are  robot  mampulators  with 
failed  actuators;  free-floating  space  robots,  where  the  base  can  be  considered  as  a  virtual 
passive  linkage  in  inertia  space;  legged  robots  with  passive  joints;  hyper-redundant  (snake¬ 
like)  robots  with  passive  joints,  etc. 

From  the  examples  above,  it  is  possible  to  justify  the  importance  of  the  study  of 
underactuated  systems.  For  example,  if  some  actuators  of  a  conventional  manipulator  fail,  the 
loss  of  one  or  more  degrees  of  freedom  may  compromise  an  entire  operation.  In  free-floating 
space  systems,  the  base  (satellite)  can  be  considered  as  a  6-DOP  device  without  positioning 
actuators.  Finally,  manipulators  with  passive  joints  and  hyper-redundant  robots  with  few 
actuators  are  important  from  the  viewpoint  of  energy  saving,  lightweight  design  and 
compactness. 

Most  of  the  le.sults  available  in  the  literature  for  fully-actuated  sys'ems  ate  difficult  to 
apply  to  these  new  ones,  because  of  the  complications  that  appear  in  their  dynamic 
formulation.  Manipulators  with  passive  joints  present  nonholonomic  constraints,  and  the 
conditions  for  integrability  of  these  constraints  are  too  stringent  [9].  In  gi  neral,  the  control 
system  must  cope  with  these  constraints  and  hopefully  take  advantage  of  tiiem  to  guarantee 
stability  and  performance  requirements.  Furthermore,  as  opposed  to  conventional 
manipulators,  in  this  problem  there  is  a  guarantee  that  no  smooth  control  law  can  achieve 
stability  of  the  system  to  an  equilibrium  point  [7],  [9].  Thus,  one  is  left  with  the  choice  of 
procuring  a  discontinuous  control  law  to  reach  a  desired  equilibrium  position,  or  being 
content  with  controlling  the  system  to  an  equilibrium  manifold. 

Some  recent  works  present  a  control  strategy  to  deal  with  the  joint  and  Cartesian  control 
problem  of  undcractuated  manipulators,  .^ai  a,".d  Tachi  [1]  presented  a  method  which 
required  brakes  lo  be  installed  on  each  passive  joints.  The  basic  idea  was  to  use  the  dynamic 
coupling  between  the  active  and  the  passive  joints,  in  order  to  drive  the  passive  joints  to  a 
desired  set-point.  Later,  Bergerman  and  Xu  [2]  enhanced  this  method  to  deal  with  parameter 
uncertainty  and  provide  the  system  with  a  greater  deal  of  robustness.  This  is  specially 
important  in  these  systems  because  the  Jacobian  mapping  from  Cartesian  to  joint  space 
depends  on  the  dynamic  parameters,  once  again  differently  from  conventional  robots,  where 
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the  Jacobian  depends  solely  on  kinematic  parameters.  Parameter  uncertainty  can  lead  the 
system  to  a  very  poor  or  even  unstable  performance.  Another  interesting  control  strategy 
requiring  the  use  of  brakes  was  done  by  Papadopoitlos  and  Dubowskj’  [10],  this  time  for  a 
space  manipulator. 

In  eveiy  work  mentioned  above,  the  authors  had  to  assume  that  “enough”  coupling 
existed  between  the  passive  and  the  active  joints,  so  that  the  controller  could  “transmit”  the 
forccs/torques  to  the  passive  joints  in  order  to  drive  them.  However,  no  attempts  were  made 
to  quantify  this  coupling,  or  to  identify  the  cases  when  it  is  too  small  as  to  be  practically 
unfeasible  to  control  the  passive  joints.  To  give  the  reader  an  introductory  feeling  of  the 
importance  of  the  dynamic  coupling  in  underactuated  mechanisms,  consider  a  Cartesian  3- 
DOF  manipulator,  where  the  joint  axes  are  mutually  perpendicular.  It  can  be  verified  that  the 
inertia  matrix  for  this  manipulator  is  diagonal  and  constant,  corresponding  to  the  physical  fact 
that  there  is  no  coupling  at  all  between  the  joints.  No  matter  how  much  one  joint  travels,  the 
other  ones  are  unaffected.  Consequently,  the  absence  of  coupling  does  not  allow  this 
mechanism  to  be  controlled  at  all  if  passive  joints  are  present. 

In  this  work,  the  authors  aim  to  provide  a  measure  of  the  dynamic  coupling  present  in 
undcractuated  systems.  Hiis  measure  is  useful  not  only  for  the  design  of  an  underactualed 
manipulator,  so  as  to  maximize  the  coupling  and  hopefully  minimize  the  energy  necessary  to 
perform  control;  it  is  also  useful  for  such  important  issues  as  actuator  placement  and  control 
strategies.  In  cases  for  which  the  number  of  actuators  is  greater  than  the  number  of  passive 
joints,  such  a  measure  can  also  be  used  in  connection  with  a  redundant  control  scheme  [8], 
in  order  to  maintain  the  system  as  far  as  possible  from  the  positions  that  yield  low  dynamic 
coupling. 


2  Dynamic  Coupling 

As  nentioned  before,  the  nonholonoraic  constraints  present  in  the  dynamic  equation  of  a 
manipulator  with  passive  joints  cannot  be  integra'.ed  in  general.  Even  partial  integrability  of 
the  acceleration  relationships  to  velocity  ones  is  not  possible  in  most  cases  [9].  This 
restriction  makes  it  impossible  to  obtain  a  direct  relationship  between  the  angles  of  the 
passive  joints  and  that  of  the  active  ones.  Thus,  it  is  nccessaty  to  work  with  the  dynamic 
equations  in  their  original  form,  and  to  try  to  derive  acceleration  relationships  to  quantify  the 
dynamic  coupling. 
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In  order  to  derive  the  measure  of  dynamic  coupling  between  the  accelerations  of  the 
passive  and  active  joints  of  an  underactuated  manipulator,  we  must  first  present  the  dynamic 
equations  governing  the  behavior  of  the  system.  Let  n  be  (he  total  number  of  joints,  r  the 
number  of  active  joints,  and  p  =  «-r  the  number  of  passive  ones.  By  using  either  the 
Newton-Euler  or  the  Lagrangian  formulation  [3],  one  can  obtain  the  following  set  of 
diffeiential  equations  relating  the  accelerations  of  the  joints  to  the  torques  supplied  by  the 
actuators; 

A/(q)4f  +  h(q,9)  =T  (1) 

Here,  the  matrix  U  is  the  it\n  inertia  matrix  of  the  manipulator,  h  is  a  vector  containing  all 
the  centrifugal,  Coriolis  and  gravitational  torques,  and  x  is  the  vector  of  torques  applied  at 
the  active  joints.  Note  that  T  has  always  p  components  equal  to  zero,  corresponding  to  the 
absence  of  actuators  at  the  passive  joints. 

Equation  (1)  is  not  useful  in  its  current  form,  for  it  does  not  reveal  the  relationship 
between  the  active  and  the  passive  joints’  velocities  and  accelerations.  If  we  paitition  the 
joint  vector  q  as: 


where  correspond  to  the  active  joint  angles  and  q^  to  the  passive  ones,  the  following 
partition  can  be  performed  on  the  dynamic  equation; 

^  ^aa  ^ap  ^ 

r  P 

It  must  be  noted  that  M  as  defined  in  (3)  i.i  not  always  equal  to  the  conventional  inertia 
matrix  of  mechanical  manipulators.  Nonetheless,  M  still  preserves  important  properties  of  the 
original  inertia  matrix,  such  as  symmetry  and  positive-definiteness.  To  see  this,  note  that  the 
new  inertia  matrix  is  obtained  from  the  original  one  after  the  swapping  of  rows  and  columns 
in  an  orderly  fashion;  if  rows  i  and  j  are  swapped,  so  must  be  columns  i  and  j. 

This  swapping  operation  can  be  represented  as  a  matrix  product.  In  order  to  avoid 
confusion,  let’s  denote  by  the  original  manipulator’s  inertia  matrix,  and  by  M  the  one 
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representing  the  active  and  passive  joints  of  the  system,  as  in  (3).  The  process  of  obtaining  M 
via  the  swapping  of  rows  and  columns  of  can  be  described  mathematically  as: 

U  =  TMJ  (4) 

where  Fisa  transformation  matrix  obtained  from  the  identity  matrix  by  the  swapping  of  rows 
i  and  j  (or  columns  i  and  j"): 


Since  T  is  obtained  from  the  identity  matrix  via  an  elementary  operation,  it  is  invertible 
(actually,  Fis  also  an  elementary  mahix).  Furthermore,  it  can  be  verified  that  it  is  equal  to  its 
inverse; 

T  =  r'*  (6) 

This  allows  us  to  write; 

M  =  TMJ'^  (7) 

and  to  establish  that  M  and  are  similar  matrices.  Now  it  is  a  known  fact  that  the  spectrum 
of  similar  matrices  are  the  same  (for  a  proof,  see  [5],  p.  152),  and  so  we  can  conclude  that  the 
new  inertia  matrix  of  the  manipulator,  M,  is  positive  definite. 

Additionally,  we  can  show  that  M  is  also  symmetric.  It  is  known  that  the  original  inertia 
matrix  is  symmetric,  and  that  the  transformation  7  is  also  symmetric.  Lhis  allows  as  to  write: 

=  TMJ  =  M  (8) 

and  to  conclude  on  the  symmetry  of  M, 
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The  submatrices  of  W  in  (3)  receive  their  indexes  according  to  the  variabies  they  relate. 
For  example,  relates  the  (null)  torques  at  the  passive  joints  to  the  acceleration  of  the 
active  ones.  The  same  reasoning  is  true  for  the  other  three  submatrices.  From  the  second  line 
of  (3),  we  can  write: 

M  q  +Af„„^  +h„  =  0  (9) 

pa^a  pp^p  p  '  ' 

or,  in  the  cases  where  Af„.  is  invertible: 

PP 

(10) 

’P  PP  po^(t  PP  P 

The  second  term  on  the  right-band  side  of  (10)  is  a  function  only  of  q  and  q ,  and  as  such  is 
completely  determined  once  measurements  of  these  variables  are  available.  Because  we  are 
focusing  on  the  acceleration  relationship  between  the  active  and  the  passive  joints,  we  rewrite 
equation  (10)  as: 

^P  “  (11) 

where 

=  ^P^K?P  (’2) 

The  acceleration  5^  can  be  viewed  as  a  virtual  acceleration  of  the  passive  joints,  generated 
by  the  acceleration  of  the  active  ones,  and  by  the  nonlinear  torques  due  to  velocity  effects. 
Given  a  desired  acceleration  of  the  passive  joints,  we  can  always  determine  at  every 
sampling  instant  the  desired  acceleration  for  as: 

^P.A0  (13) 

The  control  problem  reduces  to  finding  the  in  (11)  that  guarantees  that: 

=  ip.d  04) 

Equation  (11)  is  important  in  the  understanding  of  how  an  underactuated  sy.stcm  works. 
Torques  cun  only  be  applied  at  those  joints  which  contain  an  actuator,  or  the  active  joints. 
These  torques  produce  the  accelerations  q^,  which  indirectly  produce  the  accelerations  ij^ 
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at  the  passive  joints.  The  passive  joints’  acceierations  can  only  be  controlled  if  the  p  x  r 
matrix  possesses  a  structure  that  allows  the  actuators  torques  to  be  transmitted 
reasonably  “well”  (in  a  sense  to  be  defined  later)  to  the  passive  joints.  Thus,  the  study  of  this 
matrbt  is  of  fundamental  importance  for  the  design  and  control  of  underactuated 
manipulators. 

To  begin  the  analysis,  note  that  matrix  is  a  function  only  of  the  robot's  configuration 
q,  and  thus  is  completely  determined  based  on  the  readings  of  the  encoders  in  all  joints.  It 
does  no*,  depend  on  or  .  Thus,  equation  (11)  can  be  regarded  as  a  linear  system  of  p 
equations  Ax  =  b,  underconstrained  for  r  > p  and  overconstrained  for  r  < p . 

One  result  that  can  be  immediately  derived  from  the  structure  of  ( 1 1 )  is; 

Proposition  1  If  row  i,  1  ^  i  ^  p ,  in  matrix  contains  only  zeros,  then  the  t-th  passive 
joint  cannot  be  controlled  via  the  dynamic  coupling  with  the  active  joints. 

This  propositions  follows  from  the  fact  that,  if  has  a  line  of  zeros,  then  the  i-th  line 
in  equation  (10)  reduces  to: 


This  equality  indicates  that  the  acceleration  of  the  i-th  passive  joint  is  not  a  function  of  any 
of  the  active  joints’  accelerations,  and  thus  cannot  be  co.itrollcd  directly. 

Example  1  Consider  a  simple  two-link  manipulator  as  shown  in  Figure  1.  Joint !  rotates 
around  the  Z  axis,  while  joint  2  rotates  around  an  axis  perpendicular  to  the  first  joint  axis.  The 
inertia  matrix  M  for  this  system  is; 


M  = 


mj/^^sin  (62)^  +  /, +/2  0 

0  +7, 

2  C2  2 


(16) 


where  m,  is  the  mass  of  link  2, 1.  are  the  inertias  of  links  i  =  1 ,  2,  and  I,  is  the  distance 

2  I  Cj 

between  joint  2  and  the  center  of  gravity  of  link  2.  Two  cases  can  be  considered  here:  either 
joint  I  is  active  and  joint  2  is  passive,  or  vice-versa.  For  the  first  case,  we  .have: 
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M. 


pa 


=  m~l  +/, 

PP  2  Cj  2 


(17) 


and  therefore: 


Af^  =  0  (18) 

Equation  (17)  indicates  that  it  is  not  possible  tocontrol  via  its  coupling  with  .Thus, 
this  underactuated  system  would  not  be  useful  for  practical  purposes.  In  the  case  for  which 
joint  1  is  passive,  the  result  is  the  same: 


and  therefore: 


Af^  =  0  (20) 

We  can  conclude  that  this  mechanism’s  stracture  does  not  allow  its  passive  joint  to  be 
controlled  through  the  coupling  with  the  active  one,  whether  the  active  joint  is  joint  1  or  2. 
Note  that  the  above  statement  does  not  imply  that  the  joints  do  not  have  any  coupling  at  all. 
In  fact,  the  second  term  in  the  right-hand  side  of  ( 10)  is  generally  non-zero,  and  so  the  passive 
Joint  may  be  disturbed  for  a  given  motion  of  the  active  one.  However,  the  acceleration  of  the 
passive  joint  due  to  the  coupling  is  non-controiiable.  ■ 

Example  2  Consider  now  the  2-tink  planar  manipulator  shown  in  Figure  2.  For  this  system, 
we  have: 


Considering  joint  I  active,  and  joint  2  passive,  we  have: 

+  +^2-  ‘V  =  '”24  +  ^2 


"2'*/ '2 


(21) 


(22) 

(23) 
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Note  that  for  this  mechanism  the  structure  does  not  prevent  torque  from  being  transmitted 
from  the  active  to  the  passive  joint,  as  it  was  the  case  in  example  1.,  A  numerical 
characterization  of  this  transmission  will  be  given  in  section  3  ■ 


3  Dynamic  Coupling  Measure 

As  long  as  the  sub-matrix  is  invertible,  we  can  study  the  relationship  between  the 
accelerations  of  the  passive  and  active  joints: 

According  to  our  definition,  Af^  is  ap  x  r  matrix.  We  must  study  the  various  possibilities  that 
can  arise  depending  on  whether  there  are  more  active  or  more  passive  joints  in  the 
mechanism.  The  first  consideration  that  can  be  made  regards  the  rank  of  matrix  .  It  is 
known  that  this  rank  obeys; 


rank  (Af^)  S  min  (p,  r)  (24) 

This  fact  will  be  used  in  the  sequence. 

•  Case  1;  r<p 

Although  this  may  not  be  common,  it  may  happen  that  the  number  of  actuators  is  smaller 
than  the  number  of  passive  joints  (e.g.,  when  two  actuators  of  a  3-DOF  arm  fail),  b  this  case, 

has  maximum  rank  r,  and  equation  ( 1 1 )  has  nr  most  one  solution.  However,  this  solution 
(if  it  exists)  is  not  interesting  in  practice,  because  the  accelerations  of  p  -  r  passive  joints  will 
depend  linearly  on  the  accelerations  of  the  other  r  passive  joints.  In  other  words,  r  passive 
joints  can  be  controlled  at  every  instant,  while  the  other  p  -  r  of  them  cannot.  We  can 
conclude  that  it  is  necessary  to  have  at  least  p  actuators  in  the  underactuated  mechanism  to 
be  possible  to  control  all  p  passive  joints  independently.  Note  that  this  result  was  already 
established  by  Arai  and  Tachi  [1];  however  in  their  work,  the  authors  reached  this  conclusion 
only  after  a  study  of  the  linearized  dynamic  equations  of  the  system. 

Although  it  is  not  possible  to  control  all  p  passive  joints  when  r  <p ,  we  can  reson  to  the 
least-square  solution  in  order  to  find  the  that  generates  the  “best"  (in  a  least-squa-ms 
sense)  for  all  (or  some)  passive  joints. 
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•  Case  2:  r  =  p 

In  this  case  we  can  obtain  at  most  one  solution,  which  exists  if  the  r  x  r  matrix  is 
invertible  (or,  in  other  words,  if  both  matrices  and  are  invertible).  A  case-by-case 
pre-analysis  of  can  show  whether  it  will  be  possible  to  control  9^  using  the  acruators  at 
the  active  joints. 

♦  Case 3;  r>p 

This  is  probably  the  most  common  case,  and  certainly  the  most  interesting  one.  Here,  we 
can  obtain  at  least  one  solution  for  the  problem  of  finding  the  that  will  generate  the 
desired  9^ ,  provided  the  rank  of  matrix  is  at  least  equal  to  p.  In  the  general  case,  infinite 

solutions  can  be  found.  One  can  choose  among  these  solutions  the  one  that  provides  the 
minimum  norm  of  ,  so  as  to  save  energy,  or  effectively  make  use  of  this  redundancy  to 
accomplish  tasks  such  as  obstacle  avoidance,  actuability  maximization  [6],  etc. 

In  any  of  (he  cases  above,  it  is  useful  to  define  a  measure  of  the  dynamic  coupling  at  any 
given  instant.  For  example,  when  dealing  with  case  3,  we  can  try  to  maximize  the  coupling 
via  the  use  of  the  redundancy  present  in  the  system.  Following  [6],  [  1 1  ],  it  is  natural  to  think 
of  the  singular  values  of  ,  which  quantify  its  “degree  of  invertibility"  and  thus  its  capacity 
to  “transmit"  the  torque  from  the  active  to  the  passive  joints.  Based  on  this,  let 
o,  S  Oj  S ...  2  be  the  e  =  min  (p,  r)  singular  values  of  .  Possible  measures  of  the 
dynamic  coupling  are: 


if  r<p 

|det(Mp( 

if  r  =  p 

(25) 

if  r>p 

In  any  case, 

p,=  n®,  (26) 

I  a  1 

We  call  as  above  the  coupling  index  of  the  underactuated  manipulator.  As  will  be  shown 
in  the  sequence,  the  coupling  index  can  be  used  as  a  design  tool  for  actuator  placement. 
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desirable  robot  configuration,  or  as  a  quantity  to  be  used  on  the  real-time  control  of  the 
manipulator. 

Example  3  Let’s  retake  example  2,  and  apply  the  coupling  index  concept  to  it.  We  saw  that; 


3/,  =  - 


r  nij/j/j  cosfOj) 


1+- 


(27) 


Since  is  a  scalar,  we  have; 


P.  = 


1  +  - 


2  ^2  2 


(28) 


Let’s  adopt  the  following  parameters  for  the  quantities  above;  Wj  =  IKg,  fj  =  0.3m, 

L  =  0.1.5m,/,  =  O.lKg-m^.  Then; 

^2  ' 


a  |1+O.37cos02|  (29) 

We  see  that,  for  this  manipulator,  the  matrix  is  always  invertible,  and  thus  control  of 
the  passive  joint  via  the  dynamic  coupling  is  always  possible.  Based  on  the  present  study,  one 
can  now  pre-analyxe  the  system  in  order  to  determine  whether  or  not  control  is  possible, 
before  making  any  attempt  to  control  it.  ■ 


4  Manipulator  Design 

The  coupling  index  derived  previously  can  be  used  effectively  on  the  design  of  the 
underactuated  system  as  a  mathematical  tool  that  determines  the  optimal  actuator  placement 
of  an  underactuated  manipulator.  In  the  following  we  will  use  a  series  of  examples  to 
illustrate  its  importance. 

Example  4  If  we  consider  the  same  manipulator  as  in  example  3,  but  now  with  joint  1  as 
the  passive  joint,  we  have  the  following  results*; 


I .  The  bars  over  the  matrices  were  added  so  as  to  asoid  confusion  with  the  previous  exainpte. 
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pa 


=  '"2[ 


+  /•, 


^pp  =  Wl'c,  +'"2['l  +  (®2)]  +  ^  +  ^2 


(30) 


Therefore: 


sr,  =  - 


m2 

+  ^2 

'”l^,  +  '"2[ 

+  /,+/2 

(31) 


2  . 


Substituting  the  values  m,  =  2Kg,  /j  =  0.2Kg  ■  m  in  addition  to  the  ones  previously 
adopted,  we  have: 


Pc 


0.1225  +  0.045cose2 
”0.4575  +  0.090cose2 
0.1225 +  O.O45COS0J 
0.4575 +  O.O9OCOS02 


(32) 


Figure  3  shows  how  and  vary  as  a  function  of  Oj .  From  this  figure,  we  can  infer 
that  it  is  “easier”  for  joint  1  to  drive  joint  2  than  vice-versa,  because  of  the  greater  coupling 
available  in  the  average  for  the  lower-actuated  manipulator  than  that  for  the  upper-actuated 
one.  Thus,  the  coupling  index  indicates  that,  for  the  purpose  of  maximizing  the  dynamic 
coupling,  joint  1  should  be  the  active  one,  and  joint  2  should  be  passive.  ■ 


Note  how  this  approach  differs  from  the  one  studied  by  Lee  and  Xu  [6],  where  the  authors 
defined  the  amability  index  of  underact'iated  manipulators.  The  actuability  index  measures 
the  arbitrariness  of  the  actuator's  ability  to  cause  acceleration  at  the  end-effector.  Thus,  i! 
relates  torques  in  the  active  Joints  and  accelerations  at  the  end-effector  in  Cartesian  space,, 
while  the  coupling  index  defined  in  this  work  relates  accelerations  of  the  active  joints  to 
accelerations  of  the  passive  ones.  The  conclusions  derived  in  {6]  and  the  ones  here  should 
not  be  compared,  for  the  indexes  operate  in  different  manners.  To  be  more  specific,  the 
coupling  index  indicates  how  much  acceleration  is  possible  to  be  obtained  at  the  passive 
joints  given  limited  accelerations  at  the  active  joints.  There  is  no  attempt  to  quantify  the 
accelerations  possible  to  be  obtained  at  the  end-effector.  The  actuability  index  indicates  how 
much  acceleration  can  be  obtained  at  the  end-effector  given  limited  torques  at  the  active 
joints.  It  does  not  attempt  to  quantify  the  accelerations  at  the  passive  joints. 
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In  order  to  better  provide  the  reader  with  an  understanding  of  the  difference  between 
these  two  measures,  we  will  briefly  present  the  relationship  of  the  acceleration  of  the  end- 
effector  to  the  torque  at  the  active  joint  of  the  manipulator  in  example  4,  for  both  upper-  and 
lower-actuated  cases.  In  the  sequence,  the  subscripts  u  and  /  will  be  used  to  denote  variables 
when  the  manipulator  is,  respectively,  upper-actuated  or  lower-actuated. 

The  accelerations  of  the  end-effector  in  x  and  y  directions  can  be  easily  found  to  be: 
i  +  IjC, j)  61  +  2/20,26, 02 +  /jC, 202+  {/jJ,  +  /2S,2)6,  +/2r,2e2] 

r  2  -2  .  .  n 

y  =  -[  (li*!  + 12^12)  ^1  +  2/2^,26162  +  /2ri262-(/iCi  +  /2‘^12)  61  “ 

As  before,  we  can  consider  the  following  virtual  Cartesian  accelerations,  generated  by  the 
accelerations  at  the  joints  and  the  nonlinear  effects; 


i  ^2*12^2] 

y  =  [('iC,+/2C,2)6,+/2C,262] 


C34) 


Ignoring  the  nonlinear  effects  provenient  from  centrifugal,  Coriolis  and  gravitational  torques, 
from  (1)  we  can  write  for  the  upper-actuated  mechanism: 


4/2, 1^1+ 4^2. 2®2  =  0 

Af,,  ,0,  +Afi  2^2  = 


62  =  -(''^2,|/"2,2)®1 

0,  =  (A/2_2/dMW)3 


(35) 


where  M.  j  denotes  the  (/jj  element  of  the  original  inertia  matrix.  Substituting  (35)  into  (34) 
we  get; 


/jS,  +/2S121 


"uJJ 


M. 


2.2 


det(M)  « 


= 


/,C,  +  /2Ct2| 


1  ^2,1)1  ^2,2 

I  "W2.2jJdet(M)^'‘ 


(36) 


For  the  lower-actuated  manipulator,  the  results  ars: 


Af,  ,61  +  3/,  262  =  0 
3/2,161+3/2,262  = 


=> 


0,  =  "(3/,, 2/3/,,, >02 

62  =  (3/,  ,/dct(3/))x, 


(37) 
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L  1^1.  iJ  ^  A  ‘ 


Lc 


M 


I'-iW 


u 

“1,1 


+  /jCjj] 


.-5^11 


w, 


1,1 


M,_,jJdet(M)  ' 


(38) 


Figure  4  presents  a  comparison  of  the  nom  of  the  end-effector  acce'eration, 
||a||  =  ,  which  is  a  function  of  Oj  only.  The  full  line  represents  the  upper-actuated 

mechanism,  and  the  dashed  line,  the  lower-actuated  one.  As  we  can  see,  the  actuability  of  the 
upper-actuated  mechanism  is  always  greater  than  that  of  the  lower-actuated  one.  This  result 
contrasts  with  that  represented  by  Figure  3,  where  we  can  see  that  the  coupling  index  of  the 
lower-actuated  mechanism  is  always  greater. 

Finally,  we  borrow  here  an  example  from  [6],  to  illustrate  the  comments  above.  For  the 
same  manipulator,  the  actuability  ellipsoid  is  shown  in  Figures  S  and  6,  respectively  for  the 
lower-  and  the  upper-actuated  mechanism.  As  wc  can  see,  the  actuability  index,  which  is 
proportional  to  the  volume  (length,  in  this  case)  of  the  actuability  ellipsoid,  is  greater  for 
upper-actuated  mechanisms.  We  refer  the  reader  to  reference  [6]  for  more  detailed 
information  on  the  actuability  index. 


Example  5  Let’s  consider  now  a  “richer"  example  of  a  3-DOF  planar  manipulator  with 
rotary  joints,  as  shown  in  Figure  7.  The  inertia  matrix  is  given  by: 


M  = 


Ml,  M|2  M|3 

^21  ^22  ^23! 
^31  ^32  ^33 


(39) 


u/,+; 


'*^13  '  *^31  =  'i 

.Wjj  =  Atjj  =  '"}{‘e,*'2‘cfi)*'i 
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+/- 
j  w*  i 


For  simplicity,  let’s  adopt  the  same  par 


cuiicicib  <t>  jii  CA<tiiipic  H,  <iiiu  duupi  lui  iiiik  j  utc 


same  mass,  length  and  inertia  of  link  2.  Then; 


O7«75tO27(»C2'f00«)CjjtO060C3 
0  3225 + 0 1  J50t  j  ♦  0  OJOf  j  j  ♦  0  060c, 
OllOO+OOSOcjjtOOSOc, 


0  3225  «  0 1  SSOr,  7  0  a30c2,>  0060c,  0 1 100  *  0  030c„  »  0  (nOc, 
0  3225  6  0060c,  0ll00«0030c, 

0110060030c,  OllOO 


(40) 


Let's  assume  r  =  2,  i.e.,  we  have  two  actuators  to  be  placed  either  on  joints  1  and  2, 1 
and  3  or  2  and  3.  In  either  case,  Af^  is  a  1  x  2  matrix,  what  indicates  that  at  least  one  solution 
to  the  problem  of  finding  for  a  desired  9p  exists,  provided  that  both  elements  of  are 
not  equal  to  zero  at  the  same  instant  (one  such  possible  solution  was  demonstrated  in  [2]). 
Let’s  compute  fhe  coupling  index  for  each  case^;, 


•  Case  1;  Joints  1  and  2  are  active,  joint  3  is  passive 


=  -[l  6  0  2727(Cj,  +  Cj)  I  +02727c,j 

•  Case  2:  Joints  1  and  3  are  active,  joint  2  is  passive 


M.  =  -I 

‘J 


•  Case  3;  Joints  2  and  3  are  active,  joint  1  is  passive 

(b  3225  *  0 1350<-,  +  0.030r„  6  0  060c, 


W, 


(41) 


(42) 


(43) 


Figures  8,  9  and  10  show  the  value  of  ,  i  =  1, 2, 3  as  a  function  both  Oj  and  63 . 
Figure  1 1  shows  all  these  indexes  combined.  A  careful  consideration  of  these  figures  shows 
that,  for  most  values  of  the  joint  angles,  p^  is  the  greatest  index  of  all  three.  This  can  be 
verified  by  the  values  in  table  1.  As  we  see,  tn  none  of  the  cases  does  p  becomes  zero  (or 
“dangerously”  close  to  zero).  This  indicates  that  at  least  one  solution  will  always  exist,  no 
matter  which  joint  is  the  passive  one.  Also,  the  choice  of  joint  3  as  the  passive  joint  increases 
the  dynamic  coupling,  and  enhances  the  control  of  the  passive  joint  by  the  active  ones.  ■ 


2  Here,  ihe  indexes  1, 2  and  3  will  be  used  to  difl'crcriialc  each  case 
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Ibble  1;  Maximum,  minimum,  average  and  standard  deviation  values 
attained  by  ,  /  =  1, 2, 3 . 
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4 


i 

max^p^J 

avgjppj 

■“(PeJ 

1 

2.0021 

0.8576 

1.4211 

0.3033 

2 

1.4774 

0.6645 

1.0607 

0.2811 

3 

0.5104 

0.3912 

0.4556 

0.0374 

5  Sensitivity  Analysis 


As  explained  in  [S],  other  indexes  derived  from  the  singular  values  of  can  be  useful. 
One  of  them  is  the  condition  number  of  M^,  defined  as  the  ratio  of  the  greatest  to  the 
smallest  singular  values.  The  condition  number  is  very  useful  in  the  analysis  of  the  sensitivity 
of  equation  (1 1).  Even  in  the  cases  where  is  “big”  enough  to  guarantee  the  existence  of 
dynamic  coupling,  the  condition  number  can  indicate  that  the  relative  errors  between  the 
acceleration  of  the  active  and  the  passive  joints  is  also  too  big.  In  these  cases,  amplified  noise 
can  disturb  the  performance  of  the  mechanism. 

We  will  perform  here  a  brief  sensitivity  analysis  of  (11),  and  see  how  this  can  influence 
the  use  of  the  coupling  index.  It  is  known  that  the  norms  of  the  accelerations  in  (11)  obey: 


1<W<1 


m 


(44) 


If  noise  is  present  in  the  system  and  e.xhibits  itself  in  the  for.Ti  of  an  error  on  ,  then 
the  corresponding  error  on  obeys; 


(45) 


If  the  smallest  singular  value  of  is  too  small,  equation  (45)  shows  that  the  acceleration  of 
the  active  Joints  may  include  a  magnified  error  due  to  the  noise  present  la  ., 


Furthermore,  from  equations  (44)  and  (45)  we  can  conclude  that  the  ratio  of  the  relative 
errois  obeys: 
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or,  equivalently, 


IH 


l|/|9p|- 


(46) 


(47) 


As  we  can  sec,  the  ratio  of  the  relative  errors  is  never  smaller  than  the  inverse  of  the 
condition  number,  and  never  larger  than  the  condition  number.  If  is  too  big,  errors  will 
easily  “travel"  along  the  mechanism  making  the  control  scheme  more  difficult  and  the 
performance  worst.  These  ill-conditioned  matrices  can  render  the  control  scheme 
useless,  even  in  the  cases  where  die  dynamic  coupling  is  quite  big. 

Example  6  This  example  is  intended  to  demonstrate  the  different  cotiCiuSicns  that  can  be 
drawn  from  the  analysis  of  the  singular  values  of  the  mauix  M^.  Namely,  we  reexamine 
example  5  considering  now  as  a  measuring  index  the  condition  number  of  .  Since  is 
a  1  X  2  matrix,  it  has  only  one  singular  value;  consequently,  its  condition  number  is  always 
equal  to  1.0,  no  matter  where  the  actuators  are  placed. 


This  example  demonstrates  that  different  measures  based  on  the  same  matrix  can  lead  to 
different  conclusions:  although  the  condition  number  of  is  constant  for  any  positioning 
of  the  actuator,  the  coupling  index  is  greater  when  the  actuator  is  located  closer  to  the  base. 


A  6nal  remark  that  can  be  made  here  is  that  the  coupling  mdex  can  be  used  i.".  the  design 
of  the  underactuated  system  in  different  ways  than  it  was  used  in  example  5.  Namely,  one 
may  design  a  manipulator  where  some  pairs  of  active-passive  joints  have  maximum 
coupling,  for  easy  of  driving  of  the  passive  joints;  at  the  same  time,  the  designer  may  want  to 
minimize  the  coupling  between  specific  pairs  of  active-active,  active-passive  or  passive- 
passive  joints.  This  is  reasonable,  since  the  coupling  present  in  manipulators  is  configuration- 
varying  and  highly  nonlinear.  At  the  same  time  the  mechanism  has  maximum  coupling 
between  certain  joints  to  drive  the  passive  joints,  it  has  minimum  coupling  between  other 
joints  to  minimize  nonlinear  effects. 
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6  Implementation  Issues 


6.1  Configuration  Design 

After  going  through  the  analysis  in  section  4,  one  can  determine  the  best  actuator 
placement  so  as  to  maximize  the  dynamic  coupling  between  the  active  and  the  passive  joints. 
One  may  also  be  interested  in  finding  out  the  arm  configurations  that  yield  the  maximum 
dynamic  coupling  for  a  given  design  of  the  manipulator  and  its  workspace.  Mathematically 
this  corresponds  to  finding  the  joint  angles  q  that  maximize  the  coupling  index.  We  will 
illustrate  this  concept  with  an  example. 


Example  7  Given  the  underactuated  manipulator  in  Figure  2,  let's  suppose  an  actuator  is 
placed  at  joint  1.  As  we  saw,  is  this  case  is  given  by  (28): 

_  I  m^/,f^^cos(e^)| 


The  joint  angle  6^  that  maximizes  expression  (28)  can  be  found  via: 


mj/,/  sinCej) 


=  0 


(48) 


or: 


02  =  0  (49) 

(note  that  6^  =  it  corresponds  to  a  point  of  minimum  -  see  also  Figure  3). 

As  w(  can  see,  the  arm  configuration  that  yields  maximum  dynamic  coupling  is  the  one 
where  the  arm  is  fully  extended.  Note  how  this  differs  from  the  results  otherwise  obtained 
when  one  considers  the  classical  manipulability  measure,  introduced  by  Yoshi  kawa  [  1 2] .  For 
a  regular  fully-actuated  2-DOF  planar  manipulator  similar  to  the  one  treated  in  this  example, 
this  configuration  yields  the  minimum  of  the  manipulability  measure.  ■ 
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6.2  Individual  Joint  Coupling 

The  coupling  index,  in  addition  to  being  useful  on  the  design  of  the  underactuated  system, 
can  also  be  used  for  real-time  control.  In  [2],  Bergerman  and  Xu  demonstrated  the  feasibility 
of  driving  a  three-link  planar  manipulator  with  only  two  actuators.  The  passive  joint  and  one 
active  joint  (the  one  closest  to  the  base)  were  controlled  first.  After  the  passive  joint  reached 
its  set-point,  it  was  braked,  and  then  both  active  joints  were  controlled  to  their  desired  set- 
points.  The  active  joint  closer  to  the  base  was  chosen  to  be  controlled  first  on  a  rationale  of 
reduced  settling  time.  However,  it  may  be  the  case  that  this  joint  has  greater  coupling  with 
the  passive  joint  than  the  other  active  one,  so  it  should  be  used  to  drive  the  passive  joint,  and 
only  be  controlled  after  the  passive  joint  is  braked. 

Basically,  the  idea  here  is  to  define  an  individual  joint  coupling  index  relating  each 
passive  joint  to  each  active  one.  Working  with  these  pairs  of  active-passive  joints,  we  can 
measure  the  dynamic  coupling  existent  between  each  of  them.  Ultimately,  our  objective  will 
be  that  of  determining  which  active  joint  should  be  chosen  to  drive  each  passive  one. 

Recall  equation  (11): 


Up  =  WA 

If  we  split  the  expression  above  on  its  constituting  lines,  we  have; 

9p,  =  +  (50) 

In  order  to  study  the  coupling  between  the  /-th  passive  and  thc/-th  aclive  joint,  we  can  assign 
zero  values  to  all  except  the  j-th  active  joint’s  acceleration: 


'■l.j 


(51) 


This  equation  resembles  ( 1 1 )  and  so  we  can  define  the  coupling  index bet’^'cen  ike  i-tk passiv, 
and  the  J-th  active  joinr. 


(52) 


For  every  row  in  the  matrix  ,  the  element  (;,y)  with  greater  magnitude  will  indicate 
that  the  greatest  dynamic  coupling  for  the  passive  joint  i  comes  from  the  aci've  joint  j. 
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Example  8  We  will  review  the  example  presented  in  [2],  using  the  concept  of  individual 
joint  coupling  index  defined  above.  The  system  under  discussion  is  the  one  shown  in  Figure 
7,  where  joints  1  and  2  are  active  and  joint  3  is  passive  (according  to  the  results  in  example 
5).  For  this  system,  the  matrix  was  given  by  (41): 

(*^23  +  *^3)  '  +  ®  2727 Cjj 

Equation  (41)  can  also  be  read  as; 

=  -  [  1  +  0.2727  (Cj3  +  Cj)  ]  -  ( 1  +  O.2727C3]  9^^  (53) 

According  to  (52)  we  have: 


p  =|1  +0.2727(c33  +  C3)| 

(54) 

Pc„  =  |1+0.2727c3| 

Figures  12  and  13  present  the  values  of  and  p^  for  all  possible  combinations  of  63 
and  83 ,  Figure  14  present  both  figures  combined.  Table  2  presents  the  maximum,  minimum, 
average  and  .standard  deviation  values  for  each  of  them. 

Table  2;  Maximum,  minimum,  average  and  standard  deviation  values 
attained  by  p,  and  p,  , 


- 1 

i 

HBfc! 

min^p^J 

avg(^p,J 

std^p,,j 

It 

1.5454 

0.4546 

1.0000 

0.2728 

12 

1.2727 

. — -1 

0.7273 

1.0000 

0.1929 

For  this  system  it  is  difficult  to  tell  which  active  joint  has  greater  coupling  with  the 
passive  one.  Although  p^  attains  greater  values  than  p^  ,  it  also  attains  smaller  values; 
and  the  averages  of  both  Indexes  ate  the  same.  This  example  shows  that  the  coupling  index, 
although  useful,  may  not  provide  the  control  engineer  with  sufficient  information  to  decide 
on  which  active  joint  should  drive  the  passive  one.  Global  indexes,  based  on  the  integral  of 
the  coupling  index,  may  be  more  useful  in  this  case,  as  it  will  be  shown  in  the  sequence.  ■ 
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63  Global  Coupling  Index 

The  coupling  index  defined  previously  is  a  local  measure;  if  measures  the  “amount"  of 
dynamic  coupling  between  the  active  and  the  passive  joints  at  every  point  in  the  joint  space. 
The  greater  the  coupling  index  at  a  particular  configuration  of  the  manipulator,  the  easier  for 
the  active  joints  to  drive  the  passive  ones. 

However,  in  design  and  path  planning  problems,  one  is  more  interested  in  the  coupling 
index  in  a  global  sense,  so  as  to  measure  the  dynamic  coupling  existent  within  the  workspace 
of  the  manipulator.  In  this  case,  a  global  coupling  index  is  more  useful,  for  it  will  measure 
the  coupling  between  the  joints  at  all  points  in  the  joint  space. 

One  possible  way  of  defining  a  global  coupling  index  is  as  follows  [4]; 


e 


where  the  integrals  above  are  taken  over  the  entire  joint  space  0  6  5i"  of  the  manipulator. 
The  use  of  the  squared  value  of  the  coupling  index  is  inert  because  the  coupling  index  is 
always  non-negative;  and  this  choice  facilitates  the  problem  of  finding  a  closed-form  solution 
to  the  above  integral,  because  the  singular  values  of  a  matrix  ,4  are  equal  to  the  positive  square 

T 

roots  of  the  nonzero  eigenvalues  of  j4  M  . 

This  choice  of  the  global  coupling  index  will  take  into  account  not  only  the  local  coupling 
between  the  joints,  but  that  available  over  the  entire  joint  space,  as  the  next  examples  show. 

Example  9  We  re-analyze  now  example  S,  where  the  objective  was  to  determine  the  best 
actuator  placement  based  on  the  dynamic  coupling  between  the  joints.  There,  the  decision  to 
place  the  actuators  on  joints  1  and  2  was  based  on  the  maximum,  minimum  and  average 
values  of  for  the  various  possible  configurations.  Here  we  can  base  this  decision  on  a 
global  index,  which  already  takes  into  account  all  these  quantities. 

For  a  1  X  2  matrix,  the  unique  singular  value  is  computed  easily  as; 


A  =  [a,  ujj  =>  o  =  Ja]  +  al  (56) 

Consequently,  comparing  above  with  the  various  M  ,  r  =  !,  2, 3 ,  of  example  5,  we  have: 
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Wc  considered  in  the  following  cnlrulstions  that  the  jotnK  Caii  rotate  freely  around  their 
respective  axis  from  0  to  2)t .  In  cases  where  there  are  physical  joint  limits  this  can  be  taken 
into  account  in  the  calculation  of  the  integrals  in  (55).  For  cases  1 , 2,  and  3  studied  in  example 
5,  we  have  the  results  shown  in  table  3. 

Thble  3:  Global  coupling  index  p*,  r  =  1, 2, 3 . 


We  can  immediately  conclude  that  case  1  is  the  one  which  provides  greater  dynamic 
coupling  between  the  active  and  the  passive  joints  in  a  global  sense.  Note  that  this  is  the  same 
conclusion  as  the  one  drawn  in  example  5,  reached  in  a  much  simpler  way.  ■ 

Example  10  The  global  coupling  index  can  also  be  used  in  conjunction  with  the  individual 
joint  coupling  index  defined  in  section  6.2.  The  idea  is  to  use  the  later  instead  of  in 
equation  (35); 


p«  =9 


jpUe 


Calculating  the  above  global  individual  joint  coupling  index  for  the  manipulator  in 
example  8,  we  have  the  results  shown  in  table  4. 

As  the  result  shows,  the  global  coupling  between  the  first  and  second  active  joints  and  the 
passive  joint  is  exactly  the  same.  For  control  purposes,  then,  one  can  choose  either  active 
joint  to  dynamically  control  the  passive  joint. 
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Table  4:  Global  individual  joint  coupling  inttexes  pf  and  pf 

C,|  C. 


i 

Pc- 

It 

1.0372 

12 

1.0372 

Note,  again  that  this  conclusion  was  drawn  immediately,  as  opposed  to  example  8,  where 
no  conclusion  could  be  drawn.  ■ 


Example  11  The  coupling  index  can  also  be  useful  for  the  purpose  of  designing  the  links  of 

an  underactuated  mechanism.  Suppose  in  the  manipulator  of  Figure  2,  with  joint  1  active,  we 

2 

have  available  the  following  parameters;  =  IKg,  =  O.lKg'm  ,  Ij  =  Im. 
Additionally,  suppose  we  want  to  determine  so  as  to  maximize  the  global  coupling  index 
p* .  In  this  case  we  have; 


l^^cos  (Gj) 


0.1 +r 


1  +  70/^  +  1001^ 

1  +  20/^  +  100/? 

c,  c. 


(59) 


(60) 


Figure  15  pre.sents  the  global  coupling  index  as  a  function  of  /^  .  As  we  can  see,  for 
/j.  =  0.316m,  it  attains  the  maximum  value  p*  =  2,250.  This  example  shows  how  the 
global  coupling  index  can  be  used  for  design  issues  other  than  actuator  placement. 


7  Conclusion 

There  has  been  considerable  progress  recently  in  the  area  of  analysis  and  control  of 
underactuated  manipulators.  The  current  literature  in  this  area  however,  has  always  assumed 
that  sufficient  dynamic  coupling  was  available  between  the  active  and  passive  joints  for 
effective  control  of  the  manipulator.  Because  this  assumption  is  .not  always  valid,  it  is  vital 
for  design,  analysis,  and  control  of  these  systems  that  the  amount  of  available  coupling  be 
quantified. 
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In  this  work,  the  authors  have  proposed  a  measure  of  the  dynamic  coupling  between  the 
active  and  passive  joints  of  an  underactuated  manipulator.  The  coupling  index  indicates 
precisely  the  amount  of  coupling  available  for  the  purpose  of  controlling  the  passive  joints  of 
the  manipulator  through  the  application  of  torques  in  the  active  joints.  We  have  also  proposed 
a  global  coupling  index,  which  indicates  the  amount  of  coupling  available  over  the  entire 
workspace  of  the  manipulator. 

As  we  have  shown  through  a  series  of  illustrative  examples,  the  proposed  indices  can  be 
used  ei^ectively  within  the  contexts  of  mechanical  design,  configuration  optimization,  and 
real-time  control.  Moreover,  when  redundant  degrees  of  freedom  are  available  in  an 
underactuated  manipulator,  these  indices  may  be  used  as  local  measures  for  use  within  a 
redundant  control  scheme. 
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Figure  4!  End-effector  acceleration  for  the  mechanism  in  Figure  2. 


Figure  5:  Actuability  index  for  the  lower-actuated  mechanism  in  Figure  2. 


Figure  6:  Actuability  index  for  the  upper-actuated  mechanism  in  Figure  2. 
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Figure  7:  Three-link  planar  manipulator  with  rotary  joints. 


Figure  8;  Coupling  index  for  the  manipulator  in  Figure  7,  when  joint  3  is  passive. 
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Figure  13:  Individual  coupling  index  between  joints  2  and  3  of  the  manipulator  in  Figure  7. 


Figure  14:  Individual  coupling  indexes  for  Ihc  ntajiipulatcr  in  Figure  7. 


Figure  15:  Global  coupling  index  in  example  1 1 . 
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